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Abstract 

We present the implementation of electroweak Zjj and Wjj production via vector 
boson fusion with fully leptonic decays at NLO QCD in the Powheg framework. 
These processes represent an important background to Higgs searches in vector 
boson fusion, but they also can be seen as signal processes to study anomalous 
triple vector boson couplings as well as the impact of a central jet veto. Observables 
related to the third jet are sensitive to the parton shower which is used, a fact which 
is demonstrated by a comparison between Pythia, the standard angular-ordered 
Herwig++ shower and the new py-ordered Dipole Shower in Herwig++. 



1 Introduction 



With the discovery of a new boson at both Atlas [T] and CMS [2] we have taken one step 
closer to the understanding of electroweak symmetry breaking. To achieve this goal, one 
needs to measure the couplings of this new boson to known Standard Model (SM) parti- 
cles very precisely. Therefore, both signal and background processes must be understood 
in detail. In order to give reliable predictions for distributions at hadron colliders, it is 
mandatory to work at next-to-leading order (NLO) QCD and, if possible, interface this 
calculation with a parton shower to keep track of additional soft and collinear radiation. 
Up to now, there exist two prescriptions to match a NLO calculation to parton showers: 
Mc@Nlo [3] and Powheg pi. We will focus on the latter. The Powheg method was 
implemented in a fully flexible program named Powheg-Box [5] which equips the users 
with all subroutines needed to match their fixed-order NLO calculation to a parton shower. 



In case of the SM Higgs boson, there already exists such a Powheg-Box implementation 
for gluon fusion in association with zero jets [6], one and two jets [7] and vector boson 
fusion (VBF) |8|. VBF results in a very specific collider signature with one forward and 
one backward jet with a large rapidity gap between them which can be used to efficiently 
suppress background stemming from QCD-induced processes. To study this signature in 
data, one can look at a W or Z boson produced in VBF |9] and apply for instance central 
jet veto (CJV) techniques [IO] to this kind of processes. Since the cross section is higher for 
electroweak gauge boson production in VBF than for Higgs production, one can test the 
theoretical predictions for these processes before going to the real Higgs signal. Therefore, 
we implemented Wjj and Zjj production in VBF with subsequent leptonic decays in the 
Powheg-Box to give a NLO prediction which can be interfaced with parton showers. 



The fixed-order a s corrections to the cross section were already calculated in 11 and 



we find good agreement with the already existing electroweak Zjj-implementation in the 



Powheg-Box 12 . The QCD induced Zjj production is also part of the Powheg- 



Box 13 and can be used to test the efficiency of VBF cuts for background-suppression. 



One goal of this work is to gain experience in interfacing an existing NLO code at fixed 
order in a s with the Powheg-Box. The processes explained in detail in this publication 



offer enough complexity to study the compatibility of Vbfnlo 14 , a fully flexible parton 
level Monte Carlo program for NLO QCD corrected cross sections and distributions, and 
the Powheg-Box. The future plan is to make more processes implemented in Vbfnlo 
available in the Powheg-Box. Additionally, we turn our attention to the influence of 
the parton shower on the studied processes. To this end, we study the pr-ordered shower 



in Pythia 15 as well as the vetoed, angular-ordered shower in Herwig++ 16 and 



the new p^-ordered HERWlG++-Dipole Shower 17 , in the following just called DS++. 
From these predictions we can estimate the influence of truncation to an angular ordered 
shower. 

This paper is organized as follows: In Section [2] we review the details of the numerical 
calculation of all three processes, focusing on the subtleties of the matching between 
Vbfnlo and the Powheg-Box. In Section [3] we will give results of our calculation, 
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showered with Pythia, Herwig++ and DS + + . Conclusions are given in Section |4| 



2 Details of the Implementation 

To interface the parton-level calculation of Vbfnlo for Zjj and Wjj production via 
VBF with shower Monte-Carlo programs we use the publicly available Powheg-Box 
framework |5|. This package equips the user with all needed subroutines to go from a 



fixed-order NLO calculation in QCD to event files in the LesHouches format 18 which 
then can be interfaced with a truncated shower. To this end, the Powheg-Box asks for 
the following ingredients: 

• The Born squared matrix elements B for each partonic subprocess. The spin- 
correlated matrix elements are not needed here since there are no external gluons 
at tree level. 

• The Born color structure in the limit of a large number of colors. 



• The phase space for the Born process, see Section [272 

• The real emission squared matrix elements. 

• The finite part of the interference term between the Born and virtual amplitude. 

• The flavor structures of the Born and real emission subprocesses. We used tagging of 
the different fermion lines as described in j8j . This means that same flavor fermions 
on the upper and lower quark line internally get a different flavor (tag) to keep them 
distinct. These tags are only used to assign the possible radiation regions, which 
are searched for automatically within the Powheg-Box. 

The details of the implementation of these ingredients will follow below. 



The local subtraction terms needed to render the cross section finite are provided by the 



Powheg-Box in the FKS framework 19 . An automated check of all singular regions 
associated with one specific parton is performed and provides a good check for the flavor 
structures as well as the ratio between Born and real terms in the infrared (IR) region. 
Since the Powheg-Box offers the possibility to generate fixed leading-order (LO) and 
NLO distributions with user-defined cuts, a cross check with fixed-order calculations per- 
formed with Vbfnlo is possible and provides a strong check for the validation of the 
implementation. 



After these checks the Powheg-Box generates events in the LesHouches format from 
the PoWHEG-Sudakov factor which can be interfaced with any p^-ordered shower like 
Pythia and DS+ + or to a truncated angular-ordered shower. Since HERWIG+ + is an 
angular-ordered shower, one has to veto radiation harder than the real emission from the 
matrix element El. This option is implemented in Herwig+ + . However, one needs a so- 
called truncated shower to account for additional wide angle, soft radiation. This feature 
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is not present in the current Herwig++ release. To estimate the effect of this additional 
soft radiation, we compare Herwig++ to the pr-ordered DS + + in our analysis. 

2.1 Matrix elements 

The matrix elements were adopted from the Vbfnlo implementation explained in detail 
in 11 . Some sample diagrams for the Born and real emission contributions for Wjj- 



production are shown in Figures [T] and [2] 




Figure 1: LO diagrams contributing to Wjj production with subsequent leptonic decay. 
(a)-(d) show resonant graphs, non-resonant graphs like (e) and (/) were also included. 

When talking of Vjj-production (V = W or Z), we mean on the one hand the resonant 
production of the vector boson with leptonic decay, where off-shell effects are fully taken 



into account through a modified version of the complex mass scheme 20 with real sin 9w 



and a Breit-Wigner integration of the propagator over the whole phase space. On the 
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Figure 2: Real emission diagrams for W production via VBF: (a) and (b) are strongly 
suppressed in the VBF region and can be omitted, (c) and (d) show sample diagrams for 
gluon-induced amplitudes. Diagrams with final state gluons are obtained by crossing. 

other hand, we also take non-resonant production of the leptons into account, see Figure 
[I] (e), (/). For Zjj-production, we also take a 7* with subsequent leptonic decay into 
account. Fermion masses were set to zero throughout and 6-quarks in the initial state 
were neglected. Also, the Cabibbo-Kobayashi-Maskawa matrix was set to the unit matrix. 
This is no approximation to the calculation as long as the flavor of the jets is not tagged 
and quark masses are neglected. 



In Figure [2] (a) and (6), both (virtual) vector bosons are time-like and contain one vector 
boson which decays into a quark-anti-quark pair. This is a real emission contribution to 
WZ production and treated as a separate process. This type of diagrams is therefore not 
considered here. More details on the used approximations can be found in ITT] . 



One problem in the Powheg-Box implementation is the presence of one (Wjj) or 
even two (Zjj) t-channel photons already at LO. Consequently, the Born cross section 
is divergent if integrated over the whole phase space, though, of course, it is well defined 
with normal jet definition cuts. To avoid these singularities, amplitudes which contain 
a t-channel photon with virtuality Q 2 < 4 GeV 2 are suppressed with a large damping 



factor, as already used in 12 . At NLO, the Q 2 < 4 GeV cut affects the real emission 



contributions, as long as photon induced processes, which absorb the divergence via the 
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photon pdf, are not taken into account. In the following the missing — > VjjX piece 
is considered as a separate electroweak contribution to Vjj production and, since these 
contributions are quite small when typical VBF cuts are imposed, they were neglected. 



2.2 Phase space 

In our implementation, there exist three possibilities to evaluate the Born phase-space 
integral. The standard procedure maps the random numbers given by the integration 



routines to the physical momenta, adopted from 12 . Since the Born contributions are 
divergent in certain regions of the phase space, one can impose a cut on the two tagging 
jets, which can be changed by the user. The invariant dilepton mass for Zjj-production 
is required to be above 20 GeV by default, thus avoiding the 7* — > singularity at 
Q 2 = for massless leptons. 



As already used in 12 and originally described in [21], there also exists the possibility to 
use a Born-suppression factor F(§ n ) instead of the generation cut on the jets' pt described 
above. One possible choice for F(§ n ) is 



This factor vanishes whenever a singular region in the Born phase space <3> n is reached. 
The underlying Born kinematics are then generated in the Powheg-Box according to a 
modified B function, 

B supp = B(<S> n )F($ n ). (2.2) 

The parameters A PT = 10 GeV and k = 2 can be changed by the user. The resulting 
events have to be reweighted by a factor 1/F($ n ). 

To speed up the generation of Powheg events it is also possible to use unweighted events 
generated by Vbfnlo as phase-space generatoiQ The main advantages of this approach 
is that, first of all, the integration over the Born variables and the optimization of the 
grid with respect to the underlying Born kinematics can be omitted, since the unweighted 
events are already flat in the (Born) phase space. These unweighted events can therefore 
be seen as the perfect LO phase-space generator. Only the integration over the three 
variables of the real emission has to be handled by the integration routine. To use this 
option, unweighted events were generated using Vbfnlo. Each event i which survived 
the unweighting procedure was reweighted by the factor 



Mb $ 



2 



pdf(4 J) 



Ji 



the Born cross section over the respective numerical value of the squared Born matrix 
element including pdfs. This factor Jj is exactly the Jacobi factor of the Born phase 



lr This option can be used by setting the variable Phasespace to 2 in the powheg. input file. 
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space. A Monte Carlo integration over iV reweighted events then reproduces the Born 
cross section: 

N N 

N U N ti\M B U ( : ] ) P dff$« x 



= °LO- (2. 
For the numerical analysis shown in Section [3] we used this third method. 

2.3 Checks 

To check the implementation, all matrix elements were compared phase-space pointwise 
with the existing Vbfnlo subroutines. With this method, the evaluated couplings and 
modified routines were validated. Agreement from 11 to 15 digits was found. It was 
also verified that the subtraction and the real emission terms cancel in the singular limit. 
Another important check is the agreement of differential fixed-order NLO distributions. 
All tested distributions agree between the Powheg-Box and Vbfnlo implementation 
within statistical errors of at most 1 %. The validation of the use of unweighted events 
was done by comparing cross sections and distributions at fixed order and after event 
generation using the three different options to generate the phase space. Good agreement 
within the statistical errors was found. For Zjj production, we also compared our im- 



plementation to 12 using generation cuts in the phase space generator. Matrix elements 
were compared phase-space pointwise and agreement at the level of 10 relevant digits was 
found. Cross sections and distributions agree within the statistical uncertainties at NLO 
and after event generation. 



3 Numerical Results 

For our numerical analysis for the LHC with center-of-mass energy of 8 TeV we operate 



with the CT10 pdf set [22] with a s (M z ) = 0.11798 as implemented in the Lhapdf package 
[23] . For the calculation of the electroweak couplings we use the input parameters Mw = 
80.398 GeV, M z = 91.1876 GeV and the Fermi constant G F = 1.16637 • 10~ 5 GeV~\ 
From these parameters the total widths of the electroweak gauge bosons are calculated 
to be T z = 2.5084 GeV and Y w = 2.0977 GeV. The QED fine structure constant is 
olqed = 1/132.341 and the weak mixing angle is sin 2 9w = 0.2226. Partons are recombined 



into jets according to the ant i- fey algorithm 24 provided by the FASTjET-package 25 
with a default distance parameter R = 0.5. 

For the numerical analysis presented below we use the following inclusive cuts: 
We require that the two highest pt jets, called tagging jets, satisfy 

> 30 GeV. (3.1) 
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All observable jets, from the NLO calculation or the Shower, are demanded to have 



p TJ > 20 GeV, (3.2) 

as well as rapidity 

//, < l->- (3-3) 
To have well-observable leptons in the central region of the detector, they should obey 

p Ty i > 20 GeV and \y t \ < 2.5. (3.4) 

Since in Zjj-production, the process 7*jj —> l + l~jj is included as well, one is forced to 
impose a cut on the invariant mass of the leptons to avoid singularities: 

m u > 20 GeV. (3.5) 

All leptons should be well separated from each other and from the jets, assured by 

AR U > 0.1 and AR jt > 0.4, (3.6) 

where AR {j = - y^) 2 + (fa 



6,) 2 . 



Due to the color singlet exchange in the t-channel, the two tagging jets are widely sepa- 
rated in rapidity and usually lie in opposite detector hemispheres. Additionally, the decay 
products of the weak boson tend to be located in the rapidity gap between the two tagging 
jets. This special configuration can be used to suppress QCD backgrounds, which have a 
higher jet activity in the central detector. We therefore demand the typical VBF-cuts 

m jj > 600 GeV, Ayf/>4, y'f x y'f < and + 0.2 < Vl < - 0.2. 

(3.7) 

The factorization and renormalization scale is set to the produced vector boson's mass 



In the following we will discuss W + jj production with decay into the leptons of the first 
family. The main findings are the same for W~jj and Zjj production so only plots 
for the W + jj case will be shown. Since we are mostly interested in the effects of the 
three parton showers, PYTHIA, HERWIG++ and DS+ + , hadronisation and underlying 
event simulations were not taken into account. We used PYTHlA-version 6.4.25 with the 
Perugia 0-tune (Feb 2009) and HERWlG++-version 2.6.1a for the standard shower and 
for DS + + . The cross sections with the VBF-cuts mentioned above are shown in Table [U 
Figure [3] shows the invariant mass of the two tagging jets and the transverse momentum 
of the charged lepton for Pythia, Herwig++ and DS + + in comparison to the fixed 
order NLO prediction of Vbfnlo. As expected, the parton showers have no effect on 
these observables except for a slight change in the normalization due to the different total 
cross sections. This change comes from events which pass the cuts in a fixed order calcu- 
lation but migrate slightly by parton shower effects to phase space regions which are not 
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W+jj 


w-jj 


Zjj 


NLO 


(253.9 ±0.3) fb 


(134.4 ±0.2) fb 


(24.47 ±0.07) fb 


Vbfnlo 


(254.0 ±0.1) fb 


(134.6 ±0.1) fb 


(24.48 ± 0.02) fb 


Pythia 


(251.0 ±0.8) fb 


(131.7 ±0.5) fb 


(24.48 ±0.18) fb 


HERWIG+ + 


(249.8 ±0.8) fb 


(131.2 ±0.5) fb 


(24.08 ±0.18) fb 


DS++ 


(245.2 ±0.8) fb 


(128.0 ±0.5) fb 


(23.56 ±0.18) fb 



Table 1: Cross sections for electroweak Vjj production including VBF-cuts (3.1 3.7) with 
subsequent decay of the vector boson into the first lepton family. The NLO cross section 
was obtained with the new Powheg-Box implementation and matches the Vbfnlo 
prediction. Pythia and Herwig results include parton shower. 
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Figure 3: Differential cross-section of the invariant tagging jet mass (left) and the trans- 
verse momentum of the charged lepton (right) of the PowHEG-prediction compared to 
the fixed order curves of Vbfnlo (black solid line). The (red) dashed-dotted line shows 
the prediction of the Powheg result showered with Pythia, the (blue) dotted line cor- 
responds to Herwig±± and the (turquoise) solid line to DS-| — K The error bars show 
the statistical error of the integration, the yellow error band in the ratio plot gives the 
statistical error on the fixed-order NLO result. 

incorporated within the cuts. Also, other observables constructed from the four-momenta 
of the tagging jets or the leptons are not affected by the parton showers, the VBF signa- 
ture of the events is therefore preserved. 

Differences occur in the differential distributions of the third hardest jet, whose matrix 
elements are only LO accurate. Figure [4] shows the pr-spectrum of the third jet and its 
location relative to the tagging jets, 

y; = y J3 -(y n +y h )/2. (3.8) 

For the plot on the left hand side the cut on the transverse momentum of the third jet was 



lowered to Pt,j 3 > 1 GeV, whereas the right plot contains the usual VBF-cuts (3.1)-(3.7) 
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Figure 4: Differential distribution of the transverse momentum of the third jet with 
Ptj 3 > 1 GeV (left) and the variable y% (3.8) (right) including all cuts, comparing the 



predictions of the three parton showers to the fixed order distributions of Vbfnlo. The 
line styles are assigned as in Figure |3j 

In the low pt region, the damping of the soft divergence due to the Sudakov factor can 
be observed for all three parton showers. Between 20 and 50 GeV, Herwig+ + predicts 
more jets than Pythia, but matches the NLO prediction, whereas Pythia and DS + + 
are in good agreement. In the tail of the distributions for hard jets with p? > 75 GeV, 
all three parton showers have lower rates than the NLO distribution. This comes from 
additional hard and/or wide angle radiation which can lead to additional jets which are 
not re-clustered in the direction of the parent parton. 



Even bigger differences occur for the differential distribution of the variable y% (3.8), 
see the right plot of Figure [i With the Ay*" 9 > 4 cut, the two tagging jets peak at 
|j/* ag | w 2.7, so 1 2/3 1 < 2.7 typically corresponds to the rapidity gap between the tagging 
jets and \y$\ > 2.7 to the third jet being positioned between the tagging jets and the 
beam axis. Pythia tends to radiate more into the rapidity gap and additionally under- 
estimates the region between the tagging jets and the beam axis, whereas Herwig+-|- 
an DS + + behave the opposite way. 



This effect gets even more pronounced if one varies the distance parameter R of the 
anti-£>r algorithm or lowers the pr-cut on the third jet, see Figure [5j The differences of 
the showers are due to the fact that Pythia tends to emit more soft partons, whereas 
Herwig+-|- and DS + + preferentially emit partons in the collinear region between the 
tagging jets and the beam axis and therefore pull the third jet in that direction. If one 
lowers the pr,j 3 -cut, Figure [5] left, Pythia fills the region between the tagging jets with 
soft partons. These can recombine into the third-hardest jet which then ends up in the 
rapidity gap. However, the y% distribution for Herwig+ + and DS++ is fairly unaf- 
fected by the lower p T j 3 -cut. The shape of the curves stay the same compared to the 
NLO prediction. The collinear region is well described by the two showers, whereas the 
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Figure 5: Left: Differential 2/3-distribution (3.8) for Ptj 3 > 10 GeV. Right: Differential 
2/3-distribution for R = 0.7. The line styles are assigned as in Figure [3j 

rate in the central region is too low. 

A similar effect can be seen when the distance parameter R of the anti-fcr algorithm is 
increased. As an example the right hand side of Figure [5] shows the y% distribution for 
R = 0.7 instead of the default R = 0.5. In Pythia, the jet activity in the rapidity gap 
rises compared to the fixed-order NLO prediction since more soft and/or collinear partons 
are clustered in the third jet. This increases the possibility of the third jet to be detected. 
In contrast, Herwig++ and DS++ tend to radiate collinearly between the tagging jets 
and the beam remnant which leads to jets with high \y£\. This collinear radiation does 
not affect the shape of y% when lowering the prj 3 -cut, it does however change by increas- 
ing R (Figure |5j right). For R = 0.7, the two Herwig+ + showers produce more jets 
in the collinear region than Vbfnlo. Since Herwig++ and DS + + predict the same 
behavior for the 2/3-variable, one can conclude that the difference between Herwig+ + 
and Pythia is not caused by wide-angle, soft radiation which is included in DS+ + . 
Therefore, truncated shower effects play a minor role. The difference between the two 
HERWlG++-showers and Pythia rather seems to depend on how the available phase 
space is filled with soft and collinear radiation. 

Additionally, the jets obtained with Pythia are broader than the HERWIG++ and 
DS + + ones, which is also responsible for the different behavior of the three showers 
in the rapidity gap. This can be seen from the differential jet shape [26] p( r) of the third 
jet, which is a measure for the jet energy flow. Following Reference |27| we define the 
differential jet shape as 



p(r) 



1 

Ar 



parton^js 



PT,parton{r - Ar/2, T + Ar/2) 



PT,j, 



(3.9) 



.73 



with r ranging between 



4? and R- ^. 



PT,partan{ r ii r 2) denotes the pt of partons in an 
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Figure 6: p(r) distribution for R 
styles are assigned as in Figure |3j 

annulus between radii r x and r 2 , i.e. r x < r = \J (<f) jz - (f> pa rton) 2 + {lJj s ~ Vpartcmf < r 2- 
The sum runs over all partons which are recombined into the thirdiet. We use Ar = 0.1 
here and the normalization assures that f Q R p{r)dr = 1. In Figure 6, the averaged p{r) is 
plotted for the third jet with distance parameter R = 0.7 for different areas of the phase 
space. To distinguish the position of the third jet, we use the variable 



2/3 



\Vh ~Vh\ 



(3.10) 



The tagging jets are localized at 



0.5, l^g | < 0.5 corresponds to the rapidity gap and 



> 0.5 to the region between the tagging jets and the beam axis. 



On the left hand side of Figure |6j the differential jet shape p(r) is plotted in the whole 
allowed phase space. Clearly, Pythia produces broader jets than Herwig++ and 
DS + + . The probability to find partons with r > 0.1 which are clustered into the third 
jet are considerably higher than for the other two parton showers. The middle plot shows 
p(r) for a third jet falling in the rapidity gap between the two tagging jets, in the right 
plot the jet falls into the collinear region between the tagging jets and the beam axis. It is 
noticeable that much of the difference between HERWIG++ and Pythia stems from jets 
in the central region. This matches the observations made before and is also correlated on 
how the available phase space for additional radiation is filled: collinear radiation leads to 
narrow jets, whereas soft, wide-angle radiation is rather uncorrelated to the parton where 
it is radiated off and can therefore broaden the jet. For Pythia, this effect is large for jets 
in the central region (Figure [6j middle), whereas the differential jet shape for jets between 
the tagging jets and the beam axis is almost the same for Pythia and DS + + . In this 
region, the jet behaves almost like in an inclusive jet sample, which is reasonably well 
described by all three showers 27 . Compared to Herwig+ + , DS+ + predicts slightly 
broader jets. This can be explained by additional soft radiation due to a low IR cut-off 
on the Sudakov factor in DS+ + . By increasing this cut-off, agreement in the differential 
jet shape and the rate of third jet can be obtained. 
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Figure 7: Jet multiplicity with the standard cuts (3.1)-(3.7) for jets in the whole allowed 



phase space (left) and in the rapidity gap between the tagging jets (right). The line styles 
are assigned as in Figure [3j 



The broadening of the third jet in the central region is also the reason why, over a large 
range of px, the pT,j 3 -curve of Pythia lies below the Herwig++ prediction (see Figure 
[4j left). Pythia radiates wide-angle partons which are not clustered into the jet and 
therefore take away part of the original pt from the parent parton. In Herwig++ many 
radiated partons get clustered along the axis of the parent parton to the third jet as well as 
additional radiation from the two tagging jets. This can be seen in the prj 3 -distribution 
for large distance parameters R > 0.5, were it exceeds even the NLO prediction. To see 
the effect of radiation coming from the tagging jets also in Pythia, the distance param- 
eters R has to be increased even more. The difference between Herwig++ and DS + + 
comes from the different normalization of the two curves. As mentioned before, DS + + 
is known to radiate more soft partons which can lead to a lower rate of the third jet once 
a minimum and a maximum rapidity threshold on the tagging jets is set. 



All this has consequences for CJV techniques, since the multiplicity of jets between the 
tagging jets is quite different for Pythia, Herwig+ + and DS + + , as shown in Figure 
[7j The first three jets, which come from the hard matrix elements, are reasonably well 
described by all three showers, whereas additional jets, which solely come from the show- 
ers, show big differences. Pythia radiates off more partons which survive the jet criteria 



(3.2, 3.3, 3.6), both in the whole allowed phase space and in the rapidity gap between 



the two tagging jets. One other important conclusion of this work is that Herwig+ + , 
as a vetoed angular-ordered shower shows the same behavior as its pr-ordered sibling 
DS + + . Therefore, at least for the processes studied here, the effect of truncation can be 
neglected. 



4 Conclusions 

We implemented Wjj and Zjj production via VBF in the Powheg-Box. The Powheg 
framework allows to interface an NLO calculation with parton showers. One basic finding 
is that, as expected, the shape of the distributions of the two tagging jets and the leptons 
are mostly independent of parton shower effects. Small changes appear in the overall cross 
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sections which arise from migration of some events to phase space regions which are not 
incorporated within the cuts. In contrast, the distribution of the third jet, which is only 
LO accurate, is sensitive to the details of the parton shower in use. Dependent on the 
cuts used, the effect on third jet distributions can easily be of the order of 30 — 40%. Since 
the standard angular-ordered HERWIG++ shower and the new pr-ordered Herwig++- 
Dipole Shower are in good agreement, we expect the effects of additional wide-angle soft 
radiation, which is missing in the vetoed angular-ordered shower, to be small. However, 
there exist sizable differences between Pythia and Herwig+ + . This is due to the fact 
that Pythia predicts broader (third) jets than Herwig+ + , especially in the central 
region of the detector. These stem from soft, wide-angle radiation. In Herwig+ + , the 
third jet tends to be located in the region outside the rapidity gap due to additional 
small-angle radiation. The difference between the two HERWlG++-showers and Pythia 
seems to be caused by the filling of the available phase space for additional radiation by 
the respective shower: Pythia tends to fill the rapidity gap between the two tagging jets 
with rather soft partons, while Herwig++ leaves the rapidity gap essentially unaltered 
and radiates additional partons preferentially in the collinear region between the tagging 
jets and the beam axis. 

These differences between the three shower predictions reflect remaining uncertainties of 
available NLO predictions. They are mostly present in the distributions of the third jet, 
since it is only LO accurate, and have to be taken into account when comparing the 
predictions to data. 
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